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Abstract 



We study energy spectra, eigenstates and quantum diffusion for one- and 
two-dimensional quasiperiodic tight-binding models. As our one-dimensional 
model system we choose the silver mean or "octonacci" chain. The two- 
dimensional labyrinth tiling, which is related to the octagonal tiling, is de- 
rived from a product of two octonacci chains. This makes it possible to treat 
rather large systems numerically. For the octonacci chain, one finds singu- 
lar continuous energy spectra and critical eigenstates which is the typical 
behaviour for one-dimensional Schrodinger operators based on substitution 
sequences. The energy spectra for the labyrinth tiling can, depending on 
the strength of the quasiperiodic modulation, be either band-like or fractal- 
like. However, the eigenstates are multifractal. The temporal spreading of a 
wavepacket is described in terms of the autocorrelation function C (t) and the 
mean square displacement d(t). In all cases, we observe power laws C(t) ~ t~ s 
and d(t) ~ t". For the octonacci chain, < 5 < 1, whereas for the labyrinth 
tiling a crossover is observed from <5 = ltoO<5<l with increasing mod- 
ulation strength. Corresponding to the multifractal eigenstates, we obtain 
anomalous diffusion with < (3 < 1 for both systems. Moreover, we find that 
the behaviour of C(t) and d(t) is independent of the shape and the location of 
the initial wavepacket. We use our results to check several relations between 
the diffusion exponent (3 and the fractal dimensions of energy spectra and 
eigenstates that were proposed in the literature. 

PACS numbers: 71.23.Ft, 05.60.+W, 71.23.An, 71.30.+h 
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I. INTRODUCTION 



Since the pioneering discoveries of quasicrystals with icosahedral,0 dodecagonal,i 
decagonal,!! and octagonalE symmetry, electronic transport phenomena arguably belong to 
the most celebrated and intriguing physical properties of these intermetallic alloys.i~@ For 
instance, the electric conductivity of icosahedral quasicrystals decreases strongly with de- 
creasing temperature and with improving structural quality of the sample, and anomalous 
transport behaviour is also observed in other quantities such as thermopower or magne- 
to conductance. Stimulated by the experimental results, a lot of effort has been spent to- 
wards a better theoretical understanding of the transport phenomena in quasicrystalline 
materials. 0HHI This is also of interest from the theoretical or mathematical point of view, 
because quasicrystals as ordered aperiodic structures are intermediate between periodically 
ordered crystals and short-range ordered amorphous solids. In particular, the anomalous 
diffusion of wave packets in quasiperiodic systems has attracted wide interest. 

Multifractal eigenstates — neither extended over the system, nor exponentially local- 
ized — exist at the metal-insulator transition of the Anderson model of localization. EH In 
tight-binding models of quasicrystals, this kind of eigenstates has also been revealed.Hill! 
Generally, the energy spectra of one-dimensional (ID) quasicrystals are singular continuous.^ 
However, in higher- dimensional cases, the energy spectra can be band-like with finite mea- 
sure, fractal-like with zero band width or a mixture of partly band-like and partly fractal-like 
character .Hll 

The diffusion properties of quasicrystals are associated with the complex eigenstates and 



energy spectra stated above.l^lc2l To describe the diffusion of a wave packet initially localized 
at some site n , one usually discusses the temporal autocorrelation function@Hll@EI] 

t 

C(t) = \ J \* no (t')\ 2 dt' (1.1) 
o 

or the mean square displacement&SEil 

d{t)= (^£K-r n fmt)A (1-2) 

where *& n (t) is the amplitude of the wavefunction at time t at the nth site which is located 
at the position r n in space. Apparently, C(t) is the time-averaged probability of a wave 
packet staying at the initial site at time t, and d(t) determines the spreading of the width 
of a wave packet. 

Generally, one finds Cit) ~ t~ s , d(t) ~ t p with < 5 < 1 and < f3 < 1 for ID 
quasiperiodic systems.0@'0Tl For higher- dimensional cases, no general results are available. 
Zhong et alM observed a transition of C(t) with the increase of the quasiperiodic modulation 
strength in simple higher-dimensional Fibonacci lattices. For small patches of the octagonal 
tiling, Passaro et a/HJ found d(t) ~ t@ with < (3 < 1 even for the case of a band-like 
spectrum. However, one of us obtained C(t) ~ t^ 1 for this case after analyzing the long- 
time behaviour .0 In fact, it is quite difficult to derive the exact long-time behaviour of C(t) 
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and d(t) from the investigation of rather small systems. Therefore, a study of a large higher- 
dimensional quasiperiodic system will be significant. In this paper, we will mainly discuss 
the diffusion properties on a 2D quasiperiodic tiling related to the octagonal quasicrystals. 
The tiling is based on the octonacci chain and thus permits us to study large systems. 

Recent investigations show that the diffusion properties are connected with the mul- 
tifractality of eigenstates and energy spectra.EJtj'EJcHl It can be rigorously proven that 
the exponent 5 ruling the decay of the autocorrelation function C(t) equals the correla- 
tion dimension D 2 of the local spectral measure associated with the initial site. 00 In ID 
quasiperiodic systems, GuarneriEil analytically deduced that (3 > Dx, where D\ is the in- 
formation dimension of the spectral measure. More recently, Ketzmerick et a/.fH argued 
that (3 is also related to the multifractal properties of the eigenstates. We shall address 
the question whether these relations exist in different quasiperiodic systems, especially in 
higher- dimensional cases. 

This paper is organized as follows. In the next section, we describe the construction 
of the labyrinth tiling and its properties that are relevant for our analysis. Afterwards, in 
Sec. [TTT], we consider a tight-binding model on the labyrinth tiling and express the eigenstates 
and eigenvalues in terms of eigenstates and eigenvalues of a tight-binding Hamiltonian on 
the octonacci chain. In Sec. [TV] we show the energy spectra and multifractal eigenstates 
for both these systems. Sec. [V] describes the diffusion properties of the octonacci chain. 



The diffusion properties of the labyrinth tiling will be emphasized in Sec. [Vl[ In Sec. |VII 
we discuss the fractal dimensions of eigenstates and eigenspectra and their relation to the 
diffusion properties. Finally, we conclude in Sec. [VIII| . 



II. THE LABYRINTH TILING 

The labyrinth tilingHcJ can be considered as a subset of the octagonal quasiperiodic 
tiling,0 and vice versaB One can build it directly from the octonacci chain. In order to 
construct the labyrinth tiling, we introduce the octonacci sequence which can be produced 
by iterating the inflation rule 

6 : lZ L LSL (2 - 1} 

on the initial word w = S. The number of letters g m in the mth iterate w m = g m (w ) 
satisfies the recursion 

9m = 2# TO _i + g m -2, go = 9i = 1- (2.2) 

The numbers of letters L and S in w m are given by f m and / m _i, respectively, which fulfill 
the same recursion relation with a different initial condition 

fm = 2/ m _l + f m -2, fo = 0, fx — 1, (2.3) 

such that g m = f m + f m -x- Their ratio in the limit sequence 

lim Jul. = i im Jul. = \ (2.4) 

fm-i m - >o ° g m -i 



m— >oo 
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is given by the silver mean A = 1 + y/2 which is a root of the quadratic equation A 2 = 2 A + 1. 
As can been seen from Eq. (|2.2p , g m is odd for all m. 

Associating with the letters S and L intervals of length 1 and A, respectively, one obtains 
a linear chain C m of N m = g m + 1 sites, which is known as octonacci or silver mean chainS 
We note that all words w m obtained from the substitution rule ( |2.1| ) are palindromes, thus 
the resulting chains are symmetric under reflection. 

The labyrinth tiling can be constructed from the Euclidean product C m x C m of two such 
chains. This product is a rectangular grid, thus its vertices can be classified into even and 
odd vertices if they can be reached from the origin by an even or odd number of steps along 
the bonds, respectively. This is completely analogous to the even and the odd sublattice 
of the square lattice. Connecting the even vertices by diagonal bonds, we obtain a finite 
approximant C m of the labyrinth tiling C. The odd vertices, when connected by diagonal 
bonds, form another labyrinth tiling L* m that is dual to C m . We note that, due to the 
palindromicity of w m , C m and C* m just differ by a 90 degree rotation.il The finite labyrinth 
tiling C m consists of atoms. An example is shown in Fig. [1]. By construction, the 

labyrinth tiling is symmetric with respect to reflection at the main diagonal. Taking this 
diagonal as the x axis and the direction orthogonal to it as the y axis, the coordinates of 
the vertices of the labyrinth tiling, labeled by k, I E Z, can be written as0 

Xk,i =u k + ui (2.5a) 
Vk,i = u fc - uj (2.5b) 

where the coordinates with even values of k + I belong to C, those with odd values of k + I 
to £*. Here, 

(2.6) 



Uk = kj v 2 + 

where [x] denotes the integer closest to x. It is easy to see that the sequence of long and 
short lengths given by u k again follows the octonacci sequence, but now the two intervals 
have lengths (A ± l)/2 which again have ratio (A + 1)/(A — 1) = A. Thus, the diagonal of 
the labyrinth Xkk — 2«fc is just a y^-scaled version of the original octonacci sequence. 



III. TIGHT-BINDING MODEL 

The energy spectra for tight-binding Hamiltonians on the labyrinth tiling were in- 
vestigated by Sire.0 For properly chosen Hamiltonians, the analysis reduces to the one- 
dimensional case, and the energy spectrum can be derived directly from those of the corre- 
sponding Hamiltonian on the octonacci chain. However, the properties of eigenstates, which 
also factorize into the product of two eigenstates of the octonacci chain, were not discussed 



in Ref. [12]. Here, we follow the same route to study the eigenvalues and eigenstates. 

Consider two identical octonacci chains in the framework of a tight-binding model with 
zero on-site potentials 

ff (1 V£ M) = **Vfc? + W&? = £ (M Vf } (3.1a) 
H^ 2 ' j) = t^flf + £ m Vgf = E^^ 2 ' j) (3.1b) 
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where superscripts (1) and (2) label the two chains and the indices i and j enumerate the 
eigenfunctions ip and eigenvalues E of the two octonacci chains. Throughout the paper, we 
employ free boundary conditions, which formally corresponds to setting ip = ipN m +i = 0. 
The hopping parameters tk and ti take values according to the octonacci sequence. We 
associate tk,i = 1 to a long bond of length A and tk,i = v to a short bond of length 1, 
respectively. 

The eigenvalues of the octonacci chain are symmetric with respect to E — 0; if ip is an 
eigenstate of H with eigenvalue E, then the state ip with amplitudes 

^ = (-l) fc ^ (3-2) 

is again an eigenstate, but has an eigenvalue —E. For E — 0, the eigenvalue equation 
reduces to the recursion 

ipk+i = -—^-^fe-i (3-3) 
tk+i 

which always yields precisely two linearly independent solutions t/j* which can be chosen to 
vanish on either even or odd sites. These have the form 

vk-x = (-irvr n ^=7, vs. = o, (3.4) 

and 

^ 2 + ,-i = o, ^ = (-ir v 2 + n ( 3 - 5 ) 

s=2 l2s 

where ip± ^ and ip2 7^ are determined, up to phases, by normalization. We note that 
one has to be careful if one employs periodic boundary conditions because these, for an odd 
length of the chain, couple the even and odd sublattices of the chain. Thus, while there are 
again two states (|3.4j) and (|3.5|) for a periodic chain of even length, only a single state at 
E = exists for an odd length of the chain. 

Multiplying the two Eqs. (|3.1a|) and ( p.lbj) , we obtain 



, 2 t2s - x 



= E^E^^f (3.6) 

where we defined 

®if = 4 M) # J) (3-7) 



as an eigenfunction on the product of the two chains with eigenvalue E^-'^E^ 2 '^ . In Eq. ( p.6|) , 
only wave function amplitudes at positions (k ± 1,1 ± 1) contribute, thus the Hamiltonian 
JJv 1 ' 2 ) corresponds to hopping along the diagonals of the product grid C m x C m . The corre- 
sponding hopping parameters are the product of two hopping parameters in the octonacci 
chain and thus take values 1, v, and v 2 for diagonals of length A + 1, \/2A + 2, and A — 1, 
respectively. 
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Thus the system in Eq. (|3.6|) naturally separates into two independent sets of equations 
with k+l even or k+l odd, respectively. In this paper, we restrict our investigation to the case 
with k + I even as the other case is completely analogous. Thus, 22~( 1,2 - ) can be interpreted 
as a tight-binding Hamiltonian with zero on-site potential defined on the labyrinth tiling 
C Clearly, the eigenvalues for the labyrinth are just products of the eigenvalues for the 
octonacci chain, and all such products appear as eigenvalues because the spectra of the two 
dual labyrinth tilings C m and C* m are identical. For the corresponding eigenfunctions on 
C, we have to construct linear combinations of the product eigenfunctions $j j ( |3.7|) which 
vanish on the vertices of the dual tiling C* m . This can be done as follows. 

Suppose tpt 1 '^ and ip^ 2 '^ are normalized eigenfunctions of the octonacci chain with eigen- 
values E^-d and E&*\ respectively. Then both = ^tyfri) and = $0-$$$*) 
( p.2[ ) are eigenfunctions of H^ 1,2 ^ with the same eigenvalue E^ 1,X >E^\ where we assume 
E$) ^ and E^ ^ 0. But from Eq. (|J) we have 

*$? = (-V k+l ® { kf (3-8) 

and thus the linear combinations 

1 Jr(v) ± = _L ± (3.9) 

are normalized eigenfunctions that vanish for odd or for even values of k + l, and thus on 
C* or C, respectively. If one or both eigenvalues E^ 1 ' 1 ' and E^ 2 ^ are zero, we can make use 
of the previously discussed eigenfunctions ip + ( p.5|) and ip~ (|3.4j) to construct the desired 
wavef unctions. If one of the eigenvalues vanishes, say, without loss of generality, E^ 1 '^ = 
and E^ 2 '^ ^ 0, the appropriate four linear combinations are 

^(^,i) ±± = -^=V (1)± (V (2,i) ± ^ (2J) ) (3-10) 
where ^M 1 ) is the wavefunction of Eqs. (|3.5|) and (|3.4j) on the first chain, and we also 



used the state tjj( 2 '^ which has an energy —E^\ Clearly, the wave functions vj/^) and 
\Er(*> 3 °) have support on C, the other two linear combinations vErC*^') and + live on 

£*. Finally, for E^ = E^ 2 ^ = 0, we have four states 

where again are supported on £, and the remaining two product states 

ij/Cw anc i ijr( l .i) on £*. j n particular, this argument proves that 22 = is a (2N m — 2)- 
fold degenerate eigenvalue for the labyrinth C m . Thus we find, as for simple tight-binding 
Hamiltonians on the Penrosei! or the octagonal Ammann-Beenker tiling, a large degeneracy 
of states in the "band" center at E = 0. However, in contrast to these well-known examples 
where the degeneracy stems from certain "eonBned" stateS that occur as a consequence 
of the local topology of the tilings, the spectral measure carried by the states at 22 = 
vanishes for the labyrinth as N m — > oo, thus it is not a finite fraction of the eigenstates that 
contributes to 22 = in this case. 

In practice, having the complete knowledge of the eigenstates for the labyrinth tiling at 
our disposal, we do not need to care too much about the precise linear combinations of states 
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derived above. Since the eigenvalues Ei, i = 1, . . . , iV of the octonacci chain are symmetric 
about zero, one can obtain the set of eigenvalues of the labyrinth tiling simply as 

{EiEj | 1 < i < f , j < i) U {E l E j | f < i < N, j < i - 1} (3.12) 

where we assume that the eigenvalues of the octonacci chain are ordered as Ei > Ej for 
% > j. The corresponding eigenvectors are most easily constructed by just restricting the 
products of eigenvectors to the sites of the labyrinth C, and re-normalizing the resulting 
eigenstate. Eq. Q3.2| ) guarantees that this procedure yields the correct results, because the 
states ip and ip just differ by an alternating sign. 



IV. ENERGY SPECTRA AND WAVEFUNCTIONS 

Following the results of Sec. |TJ, one can easily calculate the density of states (DOS) 
and the integrated density of states (IDOS). For comparison, we show the DOS and the 
IDOS for the octonacci chain and the labyrinth tiling in Fig. For the octonacci chain, the 
IDOS is a devil's staircase even for v close to 1 and the DOS is singular continuous with 
zero Lebesgue measure. By more detailed analysis, one finds a self-similar energy spectrum 
for the octonacci chain with a hierarchical gap structure as described by the gap labelling 
theoremfl In contrast, we observed a smooth IDOS without visual gaps as v approaches 1 
in the labyrinth tiling. A more detailed analysis of the IDOS and the energy spectra shows 
that in the regime 0.6 < v < 1.0 the energy spectrum contains no or a finite number of gaps; 
for v < 0.6 the spectrum is fractal-like and the IDOS is similar to a devil's staircase. Sire 
found that the spectrum is singular continuous with finite Lebesgue measure for v > 0.4,0 
which may indicate that the spectrum is a mixture of band-like and fractal-like parts in 
the regime 0.4 < v < 0.6. In Fig. 0(b) one can see a peak at the center of the spectrum 
which is due to the degenerate states at E — 0. But it differs from the localized states 
observed in the Penrose tiling!! in the sense that no leap at E = is seen in the IDOS, in 
agreement with the results of the previous section. For varying parameter v, we find three 
regions with different behaviour of the DOS of the labyrinth tiling: a maximum around 
the center, distinct shoulders located between the spectral center and edge, and a tail at 
the band edge, which is similar to the behaviour observed for a tight-binding model on the 
icosahedral Ammann-Kramer tiling.0 

In order to characterize the eigenstates, we employ a multifractal analysis, which is based 
on the standard box-counting procedure. 00 In our numerical calculations, we determine the 
singularity strength a(q) and the corresponding fractal dimension f(q) by a linear regression 
procedure, but prior to this we need to check the linearity of ^/^ln//j versus hie, where 
Hi(q,e) denotes the normalized qth moment of the box probability for boxes of linear size 
eL. A homogenously extended wave function corresponds to a(q) = f(q) = d, where d 
denotes the spatial dimension. For critical eigenstates, the fractal dimension / is a smooth 
convex function of a, and a should be limited to a (/-interval. Moreover, the generalized 
dimensions of the eigenstate ip are given by D% = [f(q) — qot(q)] /(l — q) for q ^ 1 and 
Df = f(l) = a(l). 

The singularity spectra f(ci) of eigenstates for both the octonacci chain and the labyrinth 
tiling show the typical shape expected for multifractal states, thus we refrain from showing 
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these here. For the octonacci chain, the eigenstates in the "band" center are more extended 
than those at the "band" edge. In this case, the curves f{ot) become fairly narrow as v 
approaches 1. Generally, the eigenstates show stronger multifractal characteristics with 
decreasing parameter v. In contrast to the behaviour observed for the Penrose tiling ,0 for 
the labyrinth tiling we do not find that the multifractal behaviour of eigenstates becomes 
significantly stronger when moving from energies at the edge towards the center of the 
"band" . We also calculated the scaling behaviour of the inverse participation number 

p-\E,V) = J2mr)\ 4 (4-1) 

r 

with respect to the size V = L d of the system, EH i.e., 

P-\E,V) ~ V'^ E) (4.2) 

for large V. A fractal eigenstate is characterized by < 7 < 1, whereas 7 = corresponds 
to a localized state, and 7 = 1 to an extended state. In general, the scaling exponent j(E) 
depends on the energy.il Numerically, one analyzes the scaling behaviour of P^ 1 (E,V) at 
an energy E by averaging over the eigenstates within a small energy interval E ± AE/2. 
The result for eigenvectors from the center and at the lower edge of the spectrum is shown 
in Fig. |3] which corroborates the multifractal nature of the eigenstates in both systems. The 
exponent 7, given by the slope, decreases, presumably continuously, from 7 = 1 for the 
periodic case v = 1 to 7 = for v = 0. 



V. QUANTUM DIFFUSION FOR THE OCTONACCI CHAIN 

In this short section, we briefly present our numerical results of the autocorrelation 
function C(t) and the mean square displacement d(t) for the octonacci chain. Further 
discussion and comparison with the results for the labyrinth will be given below. 

Fig. |] shows the autocorrelation function C(t) of the octonacci chain. The initial site is 
located at the center of the system. The long-time behaviour of C (t) follows C (t) ~ t~ s with 
< 5 < 1 for different v. For small v, C(t) displays strong oscillatory behaviour, which may 
result from level fluctuations. The result for d(t) is displayed in Fig. |^. Evidently, d(t) ~ t@ 
and f3 increases with increasing v, limited by (3 < 1. For a given modulation parameter v, we 
observe the relation (3 > 5 between the two exponents. Similar results have been obtained 
for ID Fibonacci chains and at the mobility edge of the Harper modelJISEj Therefore, 
in accordance with the singular continuous energy spectra and the multifractal eigenstates, 
the diffusion is usually anomalous in ID quasiperiodic systems. 



VI. QUANTUM DIFFUSION FOR THE LABYRINTH TILING 

We now switch to the more interesting case of the labyrinth tiling. In Fig. ^|, we show 
the behaviour of C(t) for the labyrinth tiling. The number of sites in our system is N 2 /2 = 
19 602 2 /2 = 192 119 202, which is much larger than other 2D quasiperiodic systems discussed 
previously such as, for instance, Fibonacci latticed and the octagonal tilingE Therefore, 
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we can utilize this system to study the long-time behaviour of C(t) more accurately than 
before. Apparently, Fig. ||] again exhibits a power law behaviour C(t) ~ t~ s . By a more 
detailed analysis, we surprisingly find a transition point at v c ~ 0.6. For v < v c the slope of 
the curves decreases with decreasing v. In the regime v > v c , the behaviour of C(t) is the 
same as for a periodic system, i.e., C{t) ~ £ _1 Jil When compared to the results of Sec. [TV], 
we see that this regime corresponds to the region where one finds band-like energy spectra. 
Since 5 equals the correlation dimension D2 of the energy spectral measure JiMa 5 = 1 is 
reasonable for the case of band-like spectra. Similar to the ID case, one still has < 5 < 1 
in the regime v < v c with fractal-like or mixed spectra for the labyrinth tiling. We expect 
that this is a general result for higher- dimensional quasiperiodic systems. Furthermore, we 
find the behaviour of C(t) is independent of the initial site, which can be observed from 
the example shown in Fig. []. Of course, as our analysis is based on numerical data for a 
finite system, we cannot possibly prove the existence of a true transition point v c , because 
we cannot rule out a rapid, but continuous change in 5 around v ~ 0.6. 

The calculation of the square displacement d(t) is numerically more expensive, thus we 
restrict ourselves to a smaller system of N 2 /2 = 578 2 /2 = 167042 sites. Nevertheless, 
this is still larger than the 2D octagonal quasicrystals studied previously^ In Fig. [7|, we 
show that the long-time behaviour is described by a power law d(t) ~ t 13 . In contrast to 
C(t), we do not find a transition point for d(t) as the parameter v is varied. As for the 
octagonal tiling^ and for the octonacci chain, < (3 < 1 for the labyrinth tiling. Therefore, 
a band spectrum does not imply ballistic diffusion in quasicrystals. It can be argued that 
the exponent (3 is associated with the correlation dimension D\ of the eigenstates. H In 
ID quasiperiodic systems, or at the metal-insulator transition in the Anderson model of 
localization, the eigenstates are multifractali'@@ and < (3 < In accordance, the 

multifractal eigenstates in 2D quasicrystals may be expected to lead to anomalous diffusion 
with < (3 < 1. Possibly, ballistic diffusion can occur in 3D quasicrystals because their 
wavefunctions are more extended.0 

So far, we assumed that the initial wave packet is a 5-function, thus we start with an 
electron that is localized at a particular site n = (k , l ) and follow the spreading of its 
wave function \Ed n °} with time. This means that, in general, all eigenstates contribute to 
the time evolution because the expansion in terms of the orthonormal basis of eigenstates 

^k,l - °k,k °l,lo ~ ^AWc) M ^ D > 

and thus the entire energy spectrum is probed. For convenience, we dropped the superscripts 
± on the wavefunctions of Eqs. ( |3.9p -( ]3~TT| ), assuming that the proper linear combinations 
are used that are supported on the labyrinth C m . Note that we do not need a complex 
conjugation in Eq. (|6.1|) because the Hamiltonian is a real symmetric and we therefore 
can choose eigenvectors that form a real orthogonal matrix. In order to check for an energy 
dependence of the diffusion, we now consider different initial wave packets ty{ na '\- E —< E +—}} 
which have a finite width and are constructed as linear combinations of eigenstates from a 
certain energy window [E — AE/2, E + AE/2]. The new normalized states can be written 
as 
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where the sum ^ is restricted to the eigenstates \Er( tJ ') with eigenvalues E^ 1 '^ E^ 2 '^ £ [E — 
AE/2,E + AE/2]B Clearly, Eq. (p]) becomes Eq. (|6.1|) if the energy interval contains 
the complete spectrum, this is nothing but the usual completeness condition of the basis of 
eigenvectors. 

We numerically checked different energy windows for the octonacci chain and the 
labyrinth tiling. Due to the high DOS around E = 0, we choose smaller intervals in the 
band center. The results in Fig. || and Fig. |9| show that the long-time behaviour of C(t) and 
d(t) hardly depends on the selection of the energy window. However, it is more complex at 
small times due to the different shapes and widths of the initial wave packets. The various 
values of d(t) at the initial time reflect the width of the initial wave packet. The smaller 
the energy interval under consideration, the wider is the initial wavepacket. In practice, in 
order to avoid that the wave packet reaches the boundary too early, the energy interval may 
not be chosen too small. 



VII. DYNAMICAL SCALING AND FRACTAL DIMENSIONS 

In ID quasiperiodic systems, it is known that the inequality j3 > D\ relates the diffusion 
behaviour and the fractal properties of the energy spectrum.^ In d dimensions, this gener- 
alizes to the inequality (3 > Di/d, thus it implies a super diffusive behaviour (3 > 1/2 in two 
dimensions if D\ > 1. In Fig. |10j, the values of the exponent (3 for the octonacci chain and 
the labyrinth tiling are shown for various values of the parameter v. In all cases, we find 
that this inequality holds. Apparently, the diffusion exponents (3 for the octonacci chain 
and the labyrinth tiling are very close, which might be due to the product structure of the 
labyrinth and its wavefunctions. 

According to a conjecture by Piechon,il (3 = D-\ of the global spectral measure for 
one-dimensional quasiperiodic models with multifractal global spectral measure. In order to 
check this relation, we calculated -D_i, but it turns out that is rather difficult to extract the 
accurate values by a linear fit due to strong oscillations in the data. However, it appears 
that the relation does not hold for general parameter values v in the octonacci chain, and 
it certainly cannot be valid for the two-dimensional system as it only involves a dimension 
that characterizes the spectral measure. 

Ketzmerick et a/.il suggested an improved inequality (3 > D 2 / T)% which is numerically 
obeyed by ID quasiperiodic models.^ As can be seen from Fig. |ll], this relation applies for 
the octonacci chain as well as for the the labyrinth tiling. However, in the two-dimensional 
case the inequality is less sharp as (3 is much larger than the ratio D 2 /D 2 , in particular for 
values of the parameter v > 0.6 where the energy spectrum is smooth and D2 ~ 1. 

For multifractal wavefunctions at the Anderson transition or at quantum Hall transitions, 
one finds D\ = dD 2 for a <i-dimensional system.i!S Above, it has been demonstrated 
that D 2 = 1 for the band-like spectra in 2D quasiperiodic systems, but the corresponding 
eigenstates are multifractal with generalized dimension D 2 < 2. Although the eigenstates 
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of 2D quasiperiodic tight-binding models are similar to the critical states at the Anderson 
transition, the equality Df = dD 2 apparently does not apply to 2D quasicrystals. 

Recently, Zhong et alxB argued that one might interpret the superdiffusive behaviour in 
aperiodic systems as a ballistic behaviour in a space of effective dimension D\ , or that this 
should at least give an upper bound on the possible values of f3. In Fig. [10], we compare the 
ratio D^/d to (3. It turns out that the values of (3 and D^/d are rather close, but that there 
seems to be a systematic deviation with (3 < D^/d for small values of v and (3 > D^/d for v 
close to 1. Therefore, at least for large values of the parameter v, it appears that this bound 
does not hold. 

Finally, we also included the values of Df/d in Fig. 10, which apparently does give an 
upper bound on the values of (3 for the models under consideration. So far, this is just an 
observation, we cannot present an argument that this should hold in general. 



VIII. CONCLUSION 

In this paper, the energy spectra, wavefunctions and quantum diffusion for the octonacci 
chain and the labyrinth tiling are studied. The labyrinth tiling is based on the octonacci 
chain, which allows us to deal with very large systems. For the octonacci chain, the energy 
spectra are singular continuous and the eigenstates are critical. The energy spectra of the 
labyrinth tiling presumably are also singular continuous, but they can be band-like (i.e., 
of finite Lebesgue measure) with zero or finite gaps, a mixture of band-like and fractal 
parts, or fractal-like upon increasing the modulation strength. However, the eigenstates are 
multifractal irrespective of the value of the modulation parameter. 

The propagation of an initial wave packet is discussed in terms of the autocorrelation 
function C(t) and the mean square displacement d(t). Numerical results show that C(t) ~ 
t~ 5 and d(t) ~ t& for the octonacci chain and the labyrinth tiling. Corresponding to the 
multifractal eigenstates, we observe < (3 < 1 for both systems. In the case of fractal-like 
or mixed energy spectra and multifractal eigenstates, we find < 5 < 1. However, for a 
band-like spectrum, C(t) ~ t" 1 as in a periodic system, which causes a qualitative change 
of behaviour in C(t) for the labyrinth tiling at a parameter value v c ~ 0.6. Similar effects 
have also been observed for Fibonacci lattices^ and for the octagonal tiling.0 

We believe that the anomalous diffusion shown in d(t) and the crossover of the auto- 
correlation C(t) will be a common phenomenon in 2D quasiperiodic systems. Of course, to 
observe the crossover in C(t) one needs a parameter that allows one to continuously move 
away from the periodic case, which is not easily at hand for the most commonly investigated 
quasiperiodic model systems such as the Penrose or the octagonal tiling. Finally, we also 
studied the influence of different initial wave packets by choosing the eigenstates from vari- 
ous energy windows. The results show that the behaviour of C(t) and d(t) does not depend 
significantly on the shape and the location of the initial wave packet. 

Comparing the values of (3 with several expressions involving the fractal dimensions 
of energy spectra and eigenstates that were proposed in the literature, we find that the 
inequality (3 > D2/D2 of Ref. |2lfholds true. However, it seems that the bound (3 < /d 
proposed recently by Zhong et alBi may be violated for parameter values v close to one, i.e., 
close to the periodic case. However, we find that the weaker condition (3 < Df/d is always 
satisfied. 
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Our present work corroborates that there are strong relations between fractal properties 
of energy spectra and wavefunctions on the one hand and the exponents describing the 
quantum diffusion on the other hand. However, it appears to be difficult to find relations 
that give quantitative agreement for one- and two-dimensional aperiodic systems. Here, a 
deeper understanding of the underlying physics is desirable. Higher-dimensional systems 
constructed as products of one-dimensional systems, such as the labyrinth tiling, provide 
useful toy examples for further investigations which can, at least, be treated numerically in 
an efficient way. 
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FIGURES 
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FIG. 1. The labyrinth tiling £4. The dashed lines represent the product C4 x C4. 
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FIG. 2. The DOS (arbitrary scale) and the IDOS (bold line) for various hopping parameters v 
for (a) the octonacci chain, and (b) the labyrinth tiling. 
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(a) (b) 




In N \n(N 2 /2) 

FIG. 3. The dependence of the inverse participation numbers P^ 1 (E, V) on the size V of (a) 
the octonacci chain, and (b) the labyrinth tiling, for E = (top) and E = E m [ n (bottom). Lines 
are least-squares fits to the data, the dotted lines correspond to P~ 1 (E, V) ~ V~ x . 
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v=0.9 
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FIG. 5. The mean square displacement d(t) for the octonacci chain with N = 1394 and initial 
site no = 697. The dotted line corresponds to ballistic motion d(t) ~ t. 
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FIG. 6. The autocorrelation function C(t) for the labyrinth tiling with v = 0.1,0.2,... ,0.9. 
The system size is N 2 /2 = 19 602 2 /2 = 192119 202, and the initial site is n = (9800,9800). For 
v = 0.3, the two dashed lines correspond to different choices of the initial site no. The dotted line 
shows C(t) ~ t~ l . 
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FIG. 7. The mean square displacement d(t) for the labyrinth tiling shows d(t) ~ t 13 with 
(3 = 0.95, 0.83, 0.66, and 0.47 for v = 0.9, 0.7, 0.5, and 0.3, respectively. The dotted line corresponds 
to ballistic motion d(t) ~ t. 
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FIG. 8. The dependence of autocorrelation function C{t) on various energy windows. The 
energy windows are arbitrarily chosen from the band center (dashed line), the band edge (dotted 
line) and between the band center and the band edge (dot-dashed line), respectively. The solid 
line is the result for the full spectrum. 
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FIG. 9. The dependence of mean square displacement d(t) on various energy windows. The 
energy windows are arbitrarily chosen from the band center (dashed line), the band edge (dotted 
line) and between the band center and the band edge (dot-dashed line), respectively. The solid 
line is the result for the full spectrum. 
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FIG. 10. The exponent compared to several quantities describing the multifractal properties 
of the energy spectra and wavefunctions. The lines are included to guide the eye; the dotted lines 
correspond to the results for the octonacci chain (d = 1), the solid lines to the labyrinth tiling 
(d = 2). 
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